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We perform Monte-Carlo simulations to analyse the structure and microscopic dynamics of 
a viscous Lennard-Jones liquid coupled to a quenched reference configuration of the same 
liquid. The coupling between the two replicas is introduced via a field e conjugate to the 
overlap Q between the two particle configurations. This allows us to study the evolution of 
various static and dynamic correlation functions across the (e, T) equilibrium phase diagram. 
As the temperature is decreased, we identify increasingly marked precursors of a first-order 
phase transition between a low-Q and a high-Q phase induced by the field e. We show in 
particular that both static and dynamic susceptibilities have a maximum at a temperature- 
dependent value of the coupling field, which defines a ‘Widom line’. We also show that, in the 
high-overlap regime, diffusion and structural relaxation are strongly decoupled because single 
particle motion mostly occurs via discrete hopping on the sites defined by the reference config¬ 
uration. These results, obtained using conventional numerical tools, provide encouraging signs 
that an equilibrium phase transition exists in coupled viscous liquids, but also demonstrate 
that important numerical challenges must be overcome to obtain more conclusive numerical 
evidence. 

Keywords: Glass transition, structure of liquids, dynamic correlation functions. 


1. Introduction 

‘Cloning’ or ‘replicating’ configurations in glassy systems has a long history, which 
parallels the use of replica calculations to deal with systems defined by interac¬ 
tions containing quenched disorder [l[, such as spin glasses or disordered ferro- 
magnets UM- The relevance of replicas in the context of the glass transition of 
supercooled liquids has emerged more recently, and it originates from the proposal 
that disordered spin glass models and supercooled liquids belong to a unique uni¬ 
versality class and undergo at low temperature a random first order transition 
(RFOT) between a liquid phase and an ideal glass phase [^. Whereas this identifi¬ 
cation holds rigorously at the mean-field level 0], much remains to be understood 
regarding the applicability of these concepts in finite dimensions, where a variety 
of fluctuation effects could drastically modify the mean-field picture . 

In recent years, it has been realized that the relevance of the RFOT construc¬ 
tion can be assessed using ‘extended’ phase diagrams, where an additional control 
parameter is introduced to probe the RFOT physics in temperature regimes that 
are more easily accessible to simulations and experiments. The general idea is that 
precursors of the putative thermodynamic glass transition develop as the temper¬ 
ature is lowered in a bulk liquid, the glass phase being ‘metastable’ with respect 
to the liquid phase 0]. Therefore, the addition of well-chosen external fields could 
potentially stabilize this phase at temperatures well above the equilibrium glass 
transition and induce equilibrium phase transitions which are direct byproducts 
of the RFOT physics. In other words, external constraints might reveal growing 
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structural correlations characterizing viscous liquids approaching the glass transi¬ 
tion in a way that is not accessible to the standard tools used for simple liquids . 
Studying constrained systems thus might be useful to reveal information about 
bulk systems, and constrained phase transitions may exist even in systems where 
no finite temperature Kauzmann transition exists [8|, the two issues being logi¬ 
cally independent. Among the possible choices for such external constraints, the 
idea of pin ning the position of a set of particles has received considerable attention 
recently (^. lid], in a variety of geometries flil] fr om finite cavities [i2-14|, to fully 
random pinning 1^21] or amorphous walls (2^. . 

In this article, we analyse the case where an external coupling to a quenched 
reference configuration is introduced via a held s conjugate to the overlap Q be¬ 
tween the two copies of the^stem. This situation has been studied analytically 


in mean-held models [g, [IT], 1244261] and hnite dimensional liquids [27|] , and it was 


also studied in computer simulations [2?h31I|. The ‘annealed’ situation where the 


coupling is between two evolving copies of the system has also received attention 
from a number of groups 24, 3C, 32-13^. For the quenched coupling of interest 
in this work, the phase digram in the (e, T) plane has been established in the 
mean-held approximation [gj. In that case, an ordinary hrst-order transition line 
emerges from the bulk glass transition temperature, which separates a low-overlap 
(uncorrelated) phase from a large-overlap (localized) phase. This hrst-order tran¬ 
sition line ends at a second-order critical point which, if present, should be in the 
same universality class as the random held Ising model [24 


2 ^. 


Our central goal is to use standard computational tools to explore the high- 
temperature region of the (e,T) phase diagram where conventional Monte-Carlo 
simulations are sufficient to achieve thermalization of the coupled system. By 
analysing several static and dynamic correlation functions, we hnd a number of 
distinctive precursors of the hrst-order transition line, which allow us to dehne a 
‘Widom line’ for our system. We also hnd that the microscopic dynamics in the 
high-overlap regime is markedly different from the one of bulk supercooled liquids, 
because single particle motion occurs mostly via spatially uncorrelated, discrete 
hopping on the sites dehned by the reference conhguration. Because this diffusion 
process is very slow, we conclude that an equilibrium exploration of the two phases 
of the model at lower temperatures is impossible using conventional means. 

Our paper is organised as follows. In Sec. [2] we briehy present the model and 
our numerical methods. In Sec. [3] we analyse the static properties of the model, 
whereas Sec. 0] presents results for dynamic correlation functions. Sec. [5] closes the 
paper with some perspectives for future work. 


2. Model and numerical methods 


We performed Monte-Carlo simulations [37[| of the 
model [ 3 ^, which is a 80 : 20 binary mixture of 
Lennard-Jones potential 


Kob-Andersen glass-forming 
particles interacting via the 


Vai3{r) = 4ea/3 



( 1 ) 


where a, 13 = {A, B} denote the particles’ species. The potential is shifted to ensure 
its continuity at the cut-off distance Tc = 2.5. In the following, all quantities will be 
expressed in reduced units, selecting eAA and aAA as units of energy and distance, 
respectively, and the Boltzmann constant is ks = 1- The interaction parameters 
are aAA = 1-0, aAB = 0.8, asB = 0.88, eAA = 1-0, eAB = 1-5 and esB = 0.5. 
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They have been chosen to prevent the mixture from crystallizing at high densi¬ 
ties and low temperatures. We have studied this particular model because it is a 
well-studied and well-characterized glass-forming system . Exploratory runs 

performed using molecular dynamics methods instead of Monte-Carlo suggest that 
the dynamical results that we report do not depend qualitatively on the specihc 
choice of a microscopic dynamics. 

We studied a system of = 1000 particles in each copy. We work in three 
spatial dimensions using periodic boundary conditions at a number density p = 
N/V = 1.2, where V is the volume of the simulation box. In our Monte-Carlo 
simulations, we perform sequential single particle displacements by drawing random 
displacements within a cube of linear size Armax centered on the particle’s position. 
The proposed displacement is accepted according to the Metropolis criterion. We 
chose Armax = 0.12, which maximizes the diffusion constant of the system in the 
unconstrained mixture at T = 1.0, which represents the onset temperature for 
the bulk Lennard-Jones system. The time unit in our Monte-Carlo simulations is 
defined as N attempts to move a particle. 

To study the coupled system, we first define the overlap Qu between two con¬ 
figurations, 1 and 2, as 


Qi 2 = ^ - |ri,i - r 2 ,i|), (2) 

where 9{x) is the Heavyside function, a = 0.3 is a coarse-graining distance for 
comparing the density profiles, and denotes the position of particle i in the 
copy a G {1,2}. 

Our numerical procedure to study the coupled glassy problem is to draw a series 
of equilibrium configurations {rij at temperature Ti, which serve as reference 
configurations. We then study a second copy of the system, {r 2 }, which evolves 
at temperature T 2 and is coupled to a given reference configuration via a field e 
conjugate to the overlap between the two configurations. Only the second copy is 
allowed to evolve in this step. The total Hamiltonian thus reads 

^{ri}({r2}) = i^Lj({r2}) - eQu, (3) 

where 

^^L.l({r}) = - (4) 

o,/3 i=l j=l 

denotes the Lennard-Jones Hamiltonian for a single copy {rj, and e > 0 biases the 
coupled system towards higher values of the overlap. In our study, we use equal 
temperatures for the quenched reference configurations and the system under study, 
Ti = T 2 = T, although other choices are possible [^, [3l[. The Hamiltonian ([3]) 
contains quenched disorder, because the positions of all particles in copy 1, {ri}, 
are held fixed. Therefore, after thermal average is performed at temperature T for 
a given realisation of the disorder, we need to perform an average over independent 
configurations of the quenched copy. We found that the disorder average plays a 
negligible role for the system size and the temperature range studied in this work. 
Therefore, we only needed to average over a small number of independent quenched 
configurations to obtain accurate static and dynamic properties of the system. A 
much more demanding disorder averaging procedure would be needed at lower 
temperatures [s^, [IH . 
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Figure 1. (a) Probability distribution of the overlap P{Q) at several e values along the isotherm T = 0.8. 

(b) Typical time series of the overlap Q{t) for the same parameters (the bulk alpha-relaxation time at 
T = 0.8 is about r 2 • 10^). The overlap fluctuations are broader and slower near e* 0.7. 


Finally we find that thermalising the constrained system defined by the Hamil¬ 
tonian Q becomes increasingly difficult when e is large and/or temperature is low. 
For a given state point defined by {e,T), we consider that a system is equilibrated 
if the two following criteria are met. First, we require that the average of the over¬ 
lap evaluated over a time window of the order of 10^, displays no systematic drift. 
Second, we require that particles move an average distance of about 3 particle di¬ 
ameters. These empirical criteria are chosen so that both single particle dynamics 
and overlap fluctuations are accurately sampled in our simulations. Remark that 
even when the overlap value is large and the system explores a limited part of the 
conhgurational part close to the reference configuration, it is important to per¬ 
form very long simulations to ensure that both static and dynamic quantities are 
probed accurately. As a consequence, we have not been able to study meaningful 
parts of the (e, T) phase diagram for temperatures below T = 0.7 (our data at 
T = 0.6 barely satisfy our criteria). Note that this temperature only represents a 
modest degree of supercooling for the unconstrained system and is far above the 
mode-coupling temperature T^ct ~ 0.435, but this already represents a challenging 
situation for the coupled system. 


3. Static properties 

We start by studying the static properties of the system defined by Hamiltonian 
m, varying e and T. Intuitively, we expect that large values of e will localize the 
particles close to the sites defined by the position of the particles in the quenched 
reference configuration in order to increase the overlap Q between the two replicas. 
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In Fig. dja), we show the probability distribution of the overlap, P{Q), along the 
representative isotherm T = 0.8 for increasing values of the coupling field e. This 
temperature is slightly below the onset temperature (T ~ 1.0) of slow dynamics 
in the unconstrained system with e = 0. As the coupling e increases, we hnd 
indeed that P{Q) becomes centered around increasingly larger values of Q. As 
discussed in Sec. [H we may expect that the system undergoes a first-order phase 
transition between a low-Q phase and a high-Q phase when temperature is low 
enough. Clearly, the data in Fig. (Ha) show that the overlap increases smoothly 
with £ at this temperature and the phase transition, if present, must occur at 
lower temperatures. 

However, we notice that the distributions are narrow for both small and large 
values of e but broaden considerably at intermediate values corresponding to in¬ 
termediate overlap values, Q « 0.5. Such a broadening of the distributions for 
intermediate values of e is consistent with the idea that the system approaches the 
critical temperature from above. 

A conhrmation of this intuition is shown in Fig. (IDb) which illustrates the typical 
time evolution of the order parameter Q{t) in the course of the Monte-Carlo simu¬ 
lations. We observe that the overlap displays slow fluctuations of large amplitude 
for intermediate values of e, whereas temporal correlations and large excursions are 
strongly suppressed when the replicas are strongly coupled in the high-Q regime, 
or nearly uncorrelated at low e. These qualitative observations suggest that the 
dynamics of the coupled system display a non-trivial variation as a function of the 
coupling. The large fluctuations of the global overlap Q represent a hrst source 
of dynamic slowing down for the coupled system, which we interpret as a form 
of critical slowing down, expected in the vicinity of a second-order critical point. 
Notice that the time series shown in Fig. [1] cover a window of about 10^ timesteps, 
whereas the bulk alpha-relaxation time of the unconstrained system at the same 
temperature is r ~ 2 • 10^, which already demonstrates that accessing the prop¬ 
erties of the coupled system in conventional numerical simulation is challenging, 
even at temperatures far above the critical point. 

In Fig. [2)^a) we show the average overlap {Q) as a function of e for all studied 
temperatures, where the moments of the distribution are defined as 

(Q-)= [\qP{Q)Q^. (5) 

Jo 

We hnd that the coupling needed to localize the system in the high-overlap regime 
decreases rapidly upon decreasing the temperature. This directly implies that the 
thermodynamic driving force which allows the system to escape from a randomly 
chosen reference conhguration also decreases with the temperature. This observa¬ 
tion is at the core of the RFOT description of the glass transition [^. 

Moreover, we observe that the crossover from the low-Q regime to the high-Q 
one becomes sharper and better defined as T decreases. This is again consistent 
with the approach to a hrst-order phase transition separating a normal, \ow-Q 
phase and a localized, high-Q phase. To identify this crossover more precisely, we 
evaluate the overlap susceptibility, 

X = N[{Q^)-{Qf]. (6) 

The results are shown in Fig. [2])b). For each studied temperature, x displays a 
well-dehned maximum at a coupling e* = e*(T), which we identify as the location 
of the crossover. The locus of £*{T) in the {£,T) plane represents the equivalent 
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Figure 2. (a) Average overlap (Q) and (b) overlap susceptibility x a function of e for various tem¬ 

peratures. The crossover from uncorrelated to localized regimes occurs at lower e when T decreases, and 
becomes sharper. 


of the Widom line (4fll | for conventional phase transitions, as it represents the 
continuation of the transition line above the critical point. It should be clear that 
our interpretations are consistent with the mean-field calculations, but that our 
numerical results do not establish that an equilibrium phase transitions actually 
exists for this system. An alternative view is that the clear crossover line that we 
detect remains a smooth crossover at any finite temperatures (l9l |. Deciding which 
is the correct scenario is beyond the scope of the present work. 

The coupling parameter e in the Hamiltonian, Eq. ([3|), allows one to bias the 
system towards higher values of the overlap Q. The small value a = 0.3 used to 
evaluate the overlap leads to the formation of ‘dimers’ formed when a particle in 
the liquid is localized close to a particle in the reference configuration. We note that 
higher order clusters (such as trimers) have an enormous energy cost in the studied 
range of temperatures and couplings, and do not form. To investigate this point 
more quantitatively, we have evaluated the radial distribution functions between 
quenched and liquid replicas at various T and e, 

No. Nf, 

“ ^ i=l j=l 

In Fig. [3Ka) we show the radial distribution function between the large particles 
in the two copies, while Fig. EKb) shows correlations between large particles in the 
liquid and small particles in the reference configuration. As e increases, a well- 
defined peak develops near r ~ 0, which quantifies the increasing localisation of 
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Figure 3. Radial distribution functions between quanched and fluid replicas QAAir) QABix)-) see 
Eq. 0. Both functions develop a clear maximum near r = 0 as £ increases, reflecting the global increase 
of the overlap via the formation of ‘dimers’. Whereas the probability to form AA dimers increases with e, 
AB dimers are more numerous near e* ^ 0.7, but are suppressed at larger coupling. 


particles near the sites defined by the reference configuration. The global overlap 
Q is in fact directly related to the area under these peaks, and the pair correlation 
functions shown in Fig. [3] are the natural outcomes of integral equation approaches 
for the problem of coupled glassy systems [13, 32]. The pair correlation functions 
develop a clear minimum near r ~ 0.6, indicating the formation of tight dimers of 
liquid and quenched particles. This phenomenon resembles the formation of clusters 
in suspensions of neutral [lH and charged ultrasoft colloidal particles. Finally, 
we note that the kink at r = a is a direct result of the Heavyside function used to 
couple the two replicas and would be absent if a smoother coupling (for instance a 
Gaussian function) had been used. 

The comparison between gAAX) and gAnX) shows that the peaks at short dis¬ 
tance in these functions behave somewhat differently. First, the overall amplitude 
of the two peaks is very different, as testified by the different vertical scales used to 
represent these two functions in Fig. [3] (the amplitude at r = 0 is about 10 times 
larger for AA than for AB pairs). This shows that particles in the liquid replica 
tend to localize to a site of the reference configuration occupied by a particle of the 
same type when the coupling is strong, although the probability to form ‘defective’ 
AB dimers is not negligible. A more careful analysis of the peaks also shows that 
the peak amplitude in gAA increases monotonically with e whereas the one in gAB 
is maximum near the crossover value e* ~ 0.7. This observation suggests that the 
large fluctuations of the overlap revealed in Fig. [1] when e ~ e* likely correspond 
to a large variety of localized configurations, where dimers of various kinds can be 
formed. When the coupling becomes stronger, the system resembles increasingly 
to the reference configuration and AB dimers are suppressed, at the expense of 
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Figure 4. (a) Self-part of the intermediate scattering function Fs{k,t) evaluated at A; = 7.4 along the 

isotherm T = 0.8 and various e coupling values, (b) Mean-squared displacement 5r'^(t) for the same state 
points. The single particle dynamics slows down monotonically by increasing e, revealing the increasing 
localisation of the particles. 


AA and BB dimers. Indeed, we find that gBB{f) essentially tracks the behaviour 
of gAAir). The conclusions are also supported by the decomposition of the global 
overlap into contributions coming from the various species. 


4. Microscopic dynamics 

Close inspection of the time series of the overlap in Fig. (Hb) suggests that the 
dynamics of the liquid replica becomes increasingly sluggish upon increasing the 
coupling. In the previous section, we (tentatively) attributed this slowing down to 
the crossing of a Widom line where overlap fluctuations are very broad. However, 
we also noticed that dimers of like-particles form with increasing coupling, which 
corresponds to an increasing localization of the liquid particles on the quenched 
sites defined by the reference configuration. In this section, we will show that this 
localisation dramatically impacts the available relaxation pathways. 

To characterize the single-particle dynamics of the liquid replica, we evaluate the 
self-part of the intermediate scattering function 
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Figure 5. Evolution of the (a) structural relaxation times Ta and of (b) the diffusion coefficients D for 
large particles as a function of e for different temperatures. The dynamics slows down dramatically when 
the crossover e*{T) (shown with arrows) is crossed. Note that the change in and D~^ is quantitatively 
different, as demonstrated in Fig.^ 


and the mean-squared displacement 




(9) 


Both time correlation functions are evaluated for the large A particles (the majority 
species) along the representative isotherm T = 0.8. The results are displayed in 
Fig. m 

The unconstrained system displays practically no caging at this temperature, in 
agreement with the results of Ref. 371] on the Monte-Carlo dynamics of the same 
model. Upon increasing e, however, the dynamics slows down significantly. The 
relaxation of Fs{k,t) at long times becomes increasingly stretched and eventually 
develops a clear plateau at intermediate times for the highest values of e. A similar 
plateau is observed in the mean-squared displacement 6r‘^{t). A comparison of the 
mean squared displacement for A and B particles (not shown here) reveals that in 
the high-Q regime, the smaller B particles are somewhat less strongly localized on 
the quenched sites and can flow more easily than the large particles. The distinct 
dynamics between large and small particles is also observed in the bulk system at 
low temperatures [s^. 

As usual, we define the structural relaxation times Tq, from the condition 
Fs{k,Ta) = 1/e and the diffusion coefficients D from the Einstein relation 
limt_).oo (5r^(r) = 6Dt. In Fig. [5] we show the e-dependence of Tq and D~^ for 
A particles. For each studied temperature, we observe a marked increase of both 
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Figure 6. Decoupling of the diffusion constant and the structural relaxation is established since the product 
Dtoc increases dramatically with e. The arrows mark the crossover e*(T). 


Ta and D ^ as the coupling approaches the crossover value e* , determined from the 
analysis of the overlap susceptibility (see Sec. El). Upon further increasing e, both 
quantities continue to increase steadily, and the single particle dynamics slows down 
dramatically. Note that all data points in Fig. El are fully equilibrated according to 
the criteria detailed in Sec. [2j 

An inspection of Fig. Elsuggests that diffusion and structural relaxation both slow 
down dramatically in the high-Q regime, but their behaviour is also strongly de¬ 
coupled, as D~^ appears to change less than Tq. This is demonstrated more clearly 
in Fig. El where we show the product for the A particles. We observe that 
Dto, remains constant until £ reaches e*, above which it starts to grow signihcantly. 
Interestingly, the deviation is visible for both species, but it is less pronounced for 
the small particles (not shown). 

Such violation of the Stokes-Einstein relation, Dtq, ~ Drj ~ cst, is expected in 


systems whose dynamics is highly heterogeneous, such as supercooled liquids m, 
although other sources of decoupling have been identified in simple fluids 441. 
The coexistence of slow and fast particles in the sample leads to a decoupling 
between Tq, and D, because the former is mostly dominated by the slow particles, 
whereas the latter is dominated by the fast ones. We remark, however, that the 
Stokes-Einstein violation at large e is more pronounced than the one observed by 
decreasing temperature in the unconstrained model. This suggests that the physical 
origin of this decoupling might be somewhat different in the constrained model than 
in the bulk. 

A first hypothesis to explain the decoupling is the existence of large-scale spatial 
correlations in the dynamics [d^. This hypothesis is reasonable, since static corre¬ 
lations are promoted by the coupling field e, as testified by the maximum in the 
static susceptibility x Fig- El 

To quantify spatially heterogeneous dynamics in a more precise manner, we eval¬ 
uate the four-point dynamic susceptibility Hi-Iil, 


Xiit) = N [{fs{k,t)'^) - 


( 10 ) 


where fs{k,t) = ' (u(i) — rj(0))] represents the instantaneous value 

of the self-intermediate scattering fnnction. In supercooled liquids, X 4 (t) is char¬ 
acterized by a single peak, located around Tq, whose height xl provides a simple 
measure of the degree of heterogeneity of the dynamics (it scales roughly as the 
dynamic correlation volume). We find that the constrained system displays similar 
peaks in X4(0) and we represent the variation of xt as a function of e in Eigure[7] 
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Figure 7. Evolution of the maximum of the four-point dynamic susceptibility X 4 a function of e. Spatial 
correlations are maximum near the static crossover e*(T) (indicated by arrows). 


for all studied temperatures. 

At high temperature, Xi is essentially constant and practically independent of e. 
When T < 1.0 the behavior changes and xt varies non-monotonically as a function 
of e, with a maximum located very close to the crossover coupling e* defined from 
the maximum of the static overlap fluctuations. This shows that spatial correla¬ 
tions in the dynamics are essentially slaved to the spatial fluctuations of the static 
overlap. In particular, we find that xt decreases rapidly as e is increased beyond 
e*(T). This implies that the collective dynamics in the localized regime is actually 
strongly suppressed and that particles move instead in a spatially uncorrelated 
manner. Therefore, the strong decoupling reported in Fig. [6] does not stem from 
the coexistence of dynamically correlated domains. 

Nevertheless, it is interesting to note that the behaviour of the static and dynamic 
susceptibilities x xa is similar only when e is large, but it differs dramatically 
when e —>■ 0 where dynamic fluctuations increase rapidly as T decreases whereas 
static ones are essentially independent of temperature for the unconstrained system 
at e = 0. The fact that static and dynamic fluctuations have distinct temperature 


temper 

dependences in the moderately supercooled regime was noted before ll|, ll5|, l22 |. 


A final notable feature of our data is observed at the lowest studied temperature 
at T = 0.6 in Fig. [71 Despite the limited dynamic range we can access, the data 
reveal an additional feature, since x\ is a decreasing function of e at small e. Such a 
decrease of Xa is reminiscent of the one observed with randomly pinned particles (^ . 

For such systems, the relevant order parameter is the concentration of pinned 
particles, which plays an analogous role to the coupling e in the present model. In 
the randomly pinned systems, however, there was no trace of a subsequent increase 
of xX as a function of the concentration of pinned particles [1^ , presumably because 
the crossover line cannot be easily approached in the time window available to 
conventional computer simulations in the case of random pinning. 

The analysis of the four-point dynamic susceptibility shows that the particle 
dynamics becomes essentially uncorrelated in the high-overlap regime. Therefore, 
the explanation of the strong decoupling behavior between diffusion and alpha- 
relaxation time should be sought in the properties of single particle motion. To 
elucidate this aspect, we evalute the self-part of the van-Hove correlation function. 


Gs{r,t) = (^S{r - |ri(t) - ri(0)|)^ . 


( 11 ) 


In Fig. [51 we show the evolution of the distribution of single particle displacements 
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Figure 8. Self-part of the van-Hove correlation function Anr^Gs{T,t) for large particles for several e at 
T = 0.8. The distributions are evaluated at times t such that Fs(k^t) = 0.2. They broaden considerably 
as e increases, and the tails develop a multi-peak structure resulting from discrete particle hopping. 


for increasing e, adjusting for each value of £ the timescale such that Fs{k, t) ~ 0.2 
(this allows us to probe displacements that are large enough to better reveal the 
structure of the van-Hove function). The distributions of displacements are essen¬ 
tially featureless and close to a Gaussian for small s, but they broaden considerably 
when £ > £*{T). This considerable broadening suggests that at any given time, the 
system is composed of fast-moving and slow-moving particles that coexist in space. 
This feature alone is sufficient to account quantitatively for a strong decoupling 
phenomenon [h^. It is interesting to note that the coexistence in space of fast and 
slow particles with increasingly heterogeneous dynamics is not associated to grow¬ 
ing spatial correlations of the dynamic relaxation; this is instead a purely local 
phenomenon. 

The data in Fig. [8] not only show that the distributions broaden with e, since 
the shape of the distributions is also changing qualitatively. The distributions are 
strongly non-Gaussian at large e, and the emergence of secondary peaks is also 
obvious. These peaks suggest that the particle dynamics evolves from a continuous 
diffusive process at small e, to a slow, intermittent, hopping process at large £. 
Physically, this means that particles jump mostly from one site to another, where 
the ‘sites’ are defined by the positions of the particles in the reference configuration. 
This interpretation is harmonious with the observation of the formation of dimers 
discussed in Sec. [3l Finally, we note that particles of both species are involved in 
the hops, although the behaviour seems less pronounced for small particles, whose 
van-Hove functions are less structured than for large particles (data not shown). 
This observation contrasts with the behavior of the unconstrained model at low 
temperatures around where the jump dynamics mostly involve the smaller B 
particles. 


5. Conclusions 

We have used conventional Monte-Carlo simulations to study a viscous liquid cou¬ 
pled to a frozen reference configuration of the same liquid at the same temperature. 
Our exploration of the high-temperature portion of the (e, T) phase diagram reveals 
a ‘Widom line’ of thermodynamic and dynamic anomalies signalling the crossover 
from an uncorrelated liquid regime to a localized one, in which particles of the sys¬ 
tem are strongly constrained to reside near the sites of the reference configuration. 
These results as compatible with the existence of a line of first-order transitions at 
lower temperature terminating at a second-order critical point. Our results show 
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that computational approaches to this problem suffer from two distinct sources 
of dynamic slowing down. First, as in any phase transition, the order parameter 
(here the global overlap) displays large fluctuations near coexistence, reflecting the 
emergence of a non-trivial free-energy profile. Second, the single-particle motion 
becomes very slow in the high-overlap region because particles must hop on the 
sites of the reference configuration in a highly constrained manner in order for the 
global overlap to remain large. Altogether, this means that at least two distinct 
strategies need to be implemented in order to study the phase diagram at lower 
temperatures to overcome both the free-energy barriers near the phase transition 
and the slowing down of the particles motion in the high-overlap regime. Work is 
in progress in this direction |3ll ]. which will hopefully allow for a direct study of 
the equilibrium phase diagram of coupled viscous liquids. 
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